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ABSTRACT 

We investigate the origin and nature of the multiple sloshing cold fronts in the core of 
Abell 496 by direct comparison between observations and dedicated hydrodynamical 
simulations. Our simulations model a minor merger with a 4xlO^^M0 subcluster 
crossing A496 from the south-west to the north-north-east, passing the cluster core in 
the south-east at a pericentre distance 100 to a few 100 kpc about 0.6 to 0.8 Gyr ago. 
The gas sloshing triggered by the merger can reproduce almost all observed features, 
e.g. the characteristic spiral-like brightness residual distribution in the cluster centre 
and its asymmetry out to 500 kpc, also the positions of and contrasts across the cold 
fronts. If the subcluster passes close (100 kpc) to the cluster core, the resulting shear 
flows are strong enough to trigger Kelvin-Helmholtz instabilities that in projection 
resemble the peculiar kinks in the cold fronts of Abell 496. Finally, we show that 
sloshing does not lead to a significant modification of the global ICM profiles but a 
mild oscillation around the initial profiles. 

Key words: galaxies: clusters: intracluster medium - galaxies: clusters: individual: 
A496 - X-rays: galaxies: clusters - methods: numerical 



1 INTRODUCTION 

During their lifetimes, galaxy clusters experience a num- 
ber of major and minor mergers. Constraining the merger 
rate and understanding the properties of mergers are impor- 
tant for cosmological studies, since cluster dynamics have 
a strong impact on cluster mass determination, which is 
widely used to set constraints to cosmological parameters 
(e.g. |Rasia et al.| |2006| |Lau et al.|[2009| [Allen et al.]|2011 ). 
The investigation of individual clusters reveals insights into 
the details of the physics of galaxy clusters and their intra- 
cluster medium (ICM). Depending on the mass ratio and 
impact parameter, mergers leave observable traces in the 
clusters, among them shocks and cold fronts (see review by 
Markevitch Vikhlininll 20071. Cold fronts (CFs) appear 



as discontinuities in X-ray images and temperature maps. 
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where the brighter and denser side is also the cooler one. 
The direction of the temperature jump distinguishes them 
from shocks, where the jump is in the opposite direction. 

CFs come in two varieties: the first are merger CFs, 
which are found in merging clusters, e.g. A3667 fVikhlinin 



et al.|2001 ), the bul let cluster IE 0657-5 6 (Markevitch et al 



2002), and A2146 (Russell et al. 2010 2011). These CFs 



are contact discontinuities between the gaseous atmospheres 
of the two merging clusters. More recently, also the shocks 
associated with mergers have been discovered in some of 
these clusters (see Mark evitcti||2010 for a review) . 

CFs of the second variety form arcs around the cool 
cores of apparently relaxed clusters and have more subtle 
temperature contrasts of a factor of around 2. |Markevitch| 



et al. ( 2001 ) proposed gas sloshing as the origin of these CFs: 



while a gas-free subcluster falls through the main galaxy 
cluster, its gravitational impact during pericentre passage 
offsets the main cluster's central ICM. Consequently, the 
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offset ICM starts falling back towards the cluster centre and 
starts sloshing inside the main potential well. Usually, the 
subcluster passes the main cluster core at some distance, it 
transfers angular momentum to the ICM, and the slosh- 
ing takes on a spiral-like appearance. The details of the 
dynamics and offset mechanisms are described in |Ascasi-| 
bar & Markevitch (2006), who performed hydrodynamical 



SPH+A^-body simulations confirming this scenario. ZuHone 



et al. (2010) presented similar simulations, but concentrated 
on the heating efficiency of the sloshing process. 



Roediger et al. (2011) (Rll hereafter) presented an ex- 



tensive study of gas sloshing in the Virgo cluster. In a one-to- 
one comparison between simulations and observations, they 
showed that the sloshing mechanism can explain all proper- 
ties of the observed CPs also quantitatively, and have iden- 
tified several new characteristic features. 

Here, we continue this approach by performing a similar 
study for the cluster Abell 496. To this end, we run a set of 
numerical simulations of minor-merger induced gas sloshing 
tailored to this cluster, produce synthetic observations and 
compare them in detail to the available observations. 

The article is structured as follows: in Sect. [2] we set 
the stage by giving a comprehensive synopsis of the observed 
sloshing signatures in A496, of course including the CPs, but 
also other features. Sect.|3]briefiy introduces the simulations 
and compares simulation results and observations in detail. 
Here, we derive the most likely merger scenario and discuss 
uncertainties and possible alternatives. In Sect. ^ we discuss 
the (im)possible identification of the responsible subcluster, 
the origin of the distortions in A496's cold fronts and fur- 
ther implications. Pinally, Sect. [5] summarises our findings. 
Por the sake of clarity, the description of the observational 
data sets and the description of the simulation method and 
cluster model are collected in Appendices [A] and [B] respec- 
tively. 



2 SETTING THE STAGE: 

SLOSHING SIGNATURES IN A496 

A496 is a nearby, X-ray bright galaxy cluster (redshift 
z ^ 0.032, Dupke et al.||20Q7l [Chilingarian et al.|[2b68„ and 
refs. therein) . Por a Hubble constant of 73 km s~ Mpc~ , the 
distance of A496 is 131 Mpc, and one arcmin corresponds to 
38 kpc. Its X-ray contours are fairly regular, but the cluster 
is somewhat elongated in N-NW to S-SE direction. Hence, 
this cluster is usually classified as a relaxed cool core clus- 
ter. Its X-ray peak coincides very well with central cD galaxy 
( Dupke et aL||2007 ). 

Apart from these global characteristics, all typical fea- 
tures of gas sloshing are found in A496. We summarise them 
m Pig. [l] and give an overview of all sloshing signatures in 
the following subsections. The comparison to the simulation 
results is discussed in Sect. |3] 



2.1 Brightness edges in the X-ray image 



The Chandra observations of Dupke & White III ( 2003 ) and 
Dupke et al. ( 2007| (D07 hereafter) revealed a series of CPs 
wrapped around the cluster core in a spiral-like fashion (top 
left panel of Pig. [l]). The major CP resides about 60 kpc 
north (N) of the cluster core and was also identified in XMM 



data dChizzardi et al.||2010[ GIO hereafter). There IS a sec- 
ondary CP about 20 kpc towards the South (S) (D07, GIO) 
although D07 split this CP into two, one towards S-SW and 
one towards SE. GIO and Tanaka et al. (2006| identified a 
further CP ^ 150 kpc S of the cluster core, which is also 
evident in the Chandra image. 



2.2 Spiral-shaped brightness excess 

The structure of the cluster becomes more evident in bright- 
ness residual maps which show the ratio of the local bright- 
ness to the average value at each radius, thus highlighting 
the deviation from circular symmetry. The second panel in 
the top row of Pig. [l] displays the residual map for the clus- 
ter centre derived from the Chandra image of DO7Q The 
CPs are the outer edges of a spiral-shaped brightness ex- 
cess structure. The secondary CP about 20 kpc south of the 
centre appears as a jump in the deficit. 

The residual map in the third panel is based on the 
background and vignetting corrected XMM-Newton image 
and extends out to 500 kpc (archival data from 2007, 2008, 
see Table Al). It reveals the same structure as the one de- 
rived by [Tanaka et al. ( 2006 ) from an earlier and shorter 
XMM observation. Here we see the brightness excess spiral 
extending out to ^ 150 kpc towards the S. Having detected 
the full extent of this central structure is helpful in con- 



straining the merger geometry. Lagana et al. (2010) find the 



same spiral-shaped morphology in their substructure map 
of A496, which highlights brightness excess w.r.t. a /3-model 
fit to the cluster. 



2.3 Large-scale and general asymmetry 

Purthermore, both, Tanaka's and our large-scale residual 
map clearly show an extended brightness excess outside of 
300 kpc N-NW of the cluster core. Rll found a very sim- 
ilar feature in their simulations for the Virgo cluster and 
identified it as a typical signature of gas sloshing. Such a 
large-scale brightness excess arises a few 100 kpc from the 
cluster centre towards the direction approximately opposite 
of the pericentre of the subcluster orbit because the cluster 
centre is displaced w.r.t. the overall cluster. Consequently, 
this feature combined with the orientation of the brightness 
excess spiral indicates a subcluster orbit roughly from SW 
towards NE with the pericentre SE of the core. 

The general asymmetry of the cluster is clearly evident 
when radial profiles towards opposite directions from the 
cluster centre are compared. As an example, we do so in 
Pig. [2] for the X-ray brightness: N is compared to S and E 
to W. The profiles were averaged over azimuthal ranges of 
30° as shown in Pig. |B4| This comparison demonstrates a 
general feature typical of gas sloshing and reported already 
by Rll: profiles from opposite sides oscillate around each 
other, switching over at the CPs. This is not only true for the 
X-ray brightness, but also for temperature and metallicity. 



^ We note that the details of the residual map depend to some 
degree on the definition of the cluster centre as demonstrated in 
Appendix |B4| However, the overall structure including the spiral- 
shaped brightness excess, the northern CP and the brightness 
excess towards the S are robust features. 
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Figure 1. Sloshing signatures in A496 - qualitative comparison of observations and simulation results. Observations are shown in the top 
row, from left to right: logarithmic X-ray image for the central region from Dupke et al~| ( |2QQ7| >, replotted in our colour scale; brightness 
residuals w.r.t. azimuthal average for the central region; brightness residuals out to 500 kpc from XMM archival data (see Table \A1\ ; 
projected temperature map out to 300 kpc, from [Lovisari et al.| ( [20TT] ), replotted in our colour scale. In the X-ray and residual images, 
we mark the cluster centre. The middle and bottom rows display the same quantities in the same colour scales for our two fiducial 
simulations. In both runs, a subcluster with 4 x IO-'^^Mq and a scale radius of 100 kpc crossed the cluster on a diagonal orbit (see black 
lines in large-scale residual maps). In the middle row, the pericentre passage took place at a distance of 100 kpc to the SE 0.75 Gyr 
ago, in the bottom row at 400 kpc 1 Gyr ago. The line-of-sight is perpendicular to the orbit. The cold fronts are the outer edges of the 
brightness excess spiral. The fiducial runs reproduce nearly all features: the spiral pattern of the cold fronts and the brightness excess, the 
orientation of the brightness excess spiral, the large-scale brightness excess towards the N-NW, the strongest large-scale brightness deficit 
towards the E, and the appearance of the temperature map. These two best simulations differ in details. The larger impact parameter 
in the bottom row leads to a more regular structure of the cold fronts, and reproduces the cold front positions towards all directions 
simultaneously (see Fig. [s] and Sect. |3.4| . The smaller impact parameter in the middle row leads to stronger perturbation in the cluster 
core and higher shear velocities along the cold fronts. While this run does not match the cold front radii towards the East and West as 
well, it resembles the disturbed shape of the observed cold fronts, i.e. there are kinks of similar size as well as a double cold front towards 
the West. 
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Figure 2. Comparison between X-ray brightness profiles from 
opposite sides. The directions are colour-coded and labelled. For 
clarity, the profiles towards E and W have been multiplied by 
0.3. Solid lines are observed profiles (XMM-Newton, GIO), dashed 
lines are simulation results for our fiducial run "distant". Dotted 
lines indicate our initial profiles. Profiles from opposite sides os- 
cillate around each other, switching over at the cold fronts. The 
vertical double lines bracket the cold fronts towards the individual 
directions, using the same colour code as the profiles. 



The large-scale asymmetry appears as the split-up between 
opposite profiles at radii outside the outermost CF. 



2.4 Temperature and metallicity structure 

The XMM-Newton temperature map (top right panel in 
Fig. [l] [Lovisari et al. 2011 Lll hereafter) confirms that 
brightness excess regions correspond to cooler regions, de- 
spite the fact that the temperature map is very coarse. 
Also the metallicity map (Lll) shows a rough correspon- 
dence of brightness excess and higher metallicity, although 
it suffers from considerable systematic uncertainties (see Ap- 
pendix [A|. The less deep Chandra observation of DOT does 
not show significant structures in the metallicity map, but 
a clear cool region just inside the northern CF in also found 
in the Chandra data. 



2.5 Cold front positions from radial profiles 

Radial profiles of X-ray brightness, temperature and metal- 
licity are shown in Fig. [3] and are discussed in compari- 
son to the simulations in Sect. 13.51 Different observational 
datasets are available and are compared in Appendix [A] 
Figure |3] shows a representative selection based on XMM- 
Newton observations: the X-ray brightness and temperature 
profiles from GIO, and our re-calculation of temperature and 
metallicity profiles with three assumptions (APEC/VAPEC 
model, hydrogen column density fixed/free during fit; see 
Appendix[A| . The different temperature and metallicity pro- 
files are shown as grey-shaded bands, together they indi- 
cate the amount of systematic error. In Fig. |Al| the different 
datasets can be distinguished. The profiles are averaged az- 
imuthally over 30^^ as shown in Fig. |B4[ 

In principle, CFs appear in radial profiles as downward 



Table 1. Positions of cold fronts in A496 derived from radial pro- 
files: we give a minimal extent in azimuth and the inner and outer 
edge of the cold front in radius in the sectors labelled in Fig. |B4| 
Also for the azimuth we follow the convention from Fig. |B4| 



azimuthal 
extent 



radius in kpc 
(in arcsec) 
in sectors 



major CF N-W: 


-10° 


to 


N: 50-63 


W: 35-60 




110° 




(80-100) 


(55-95) 


secondary CF E-S-SW: 


180° 


to 


S: 13-27 


E: ?-15 




300° 




(20-43) 


(?-25) 


outer CF E-S-SW: 


180° 


to 


S: 114-177 


E: 82-184 




300° 




(180-280) 


(130-290) 



jumps in X-ray brightness and metallicity, and as upward 
jumps in temperature. The position of these jumps identifies 
the CF position or CF radius. However, projection and the 
combination of azimuthal and radial binning reduces the 
intrinsic discontinuities to more or less steep gradients. 

In the radial profiles towards the N, the major, northern 
CF is clearly identified by such a steep gradient in all quan- 
tities. The other CFs are more subtle because they reside in- 
side the cool core near the cluster centre. Here, not only the 
overall X-ray brightness and metallicity are strongly declin- 
ing but, additionally, the overall temperature increases with 
radius. This means that the general gradients of all quanti- 
ties are in the same direction as typical for a CF. Thus, the 
gradient due to the CF blends with the overall gradient and 
the CF radius is hard to identify. As a consequence, all but 
the major northern CF are mainly apparent as kinks in the 
observed X-ray brightness profiles, and their exact position 
is difficult to constrain from the radial profiles. 

Here, the oscillating behaviour of opposite profiles as 
described in Sect. |2.3| and Fig.[2]can help: the radius at which 
the profiles switch over identifies the position of the CF. 
Combining the information from each profile alone and from 
the comparison of opposite directions, we infer the positions 
of the CFs as listed in Table [T] For each CF we state its 
azimuthal range determined from the maps in Fig. [l] and 
a radial range which encompasses the switch-over region as 
seen in Fig. [2] For the secondary CF towards the E, the 
azimuthal profiles do not allow to reliably constrain the inner 
edge of the CF. We measure the CF width more accurately 
in the next subsection. 

Our results mostly agree with the ones from D07 and 
GIO. We have extended the major (northern) CF towards 
the W and the outermost southern CF towards the E. We 
find the outer southern CF at the same radius as GIG. We 
have also interpreted the two inner southern CFs from D07 
as one secondary CF. The same CF was reported by GIO at 
a slightly larger radius and smaller azimuthal extent. 

We mark these CF positions in all relevant plots. 



2.6 Detailed structure and width of the cold 
fronts 

The inner CFs in A496 are peculiar because of their 
boxy morphology. Instead of smooth arcs as e.g. in Virgo 
dSimionescu et al [ [2010) or A2142 ( [Markevitch et al.||2000" 



Owers et al.|2009f they show several kinks. We discuss their 
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Figure 3. Comparison of observed and synthetic profiles of X-ray brightness (top panels), projected temperature (middle panels) and 
metallicity (bottom panels) towards N and S (Ihs figure), and E and W (rhs figure). For the observed X-ray profiles, we plot only every 
4th available data point for the sake of clarity. The profiles are extracted from 30° azimuthal sectors as shown in Fig. |B4| We show 
the following XMM-Newton observational data: X-ray brightness and temperature from GIO, and our re-calculation of temperature and 
metallicity profiles with different assumptions (APEC/VAPEC model, hydrogen column density fixed/free during fit; see Appendix [a}. 
The different resulting temperature and metallicity profiles are shown as grey-shaded bands, together they indicate the amount of 
systematic error (see Fig. |Al| for a distinction). We plot three representative simulations: the fiducial "distant" run, the fiducial "close" 
run, and a low mass case with a subcluster of 10^^ Mq and a pericentre distance of 100 kpc; see Table B2 for model parameters. The 



dashed vertical double lines mark the position of the cold fronts in the radial profiles (just single line for the inner eastern cold front) 



possible origin in Sect. |4.4l The CPs being non-circular, ra- 
dial azimuthally averaged profiles are not optimal for mea- 
suring their widths. Hence, we derived additional surface 
brightness profiles across the N and W cold front and the 
outer southern front, taking particular care to place the pro- 
file directions locally perpendicular to the fronts. For this 
reason, our profiles are not extracted from annuli but the 
box-like regions shown in the left-hand part of Fig. [4] The 
right-hand-side of this figure compares the resulting bright- 
ness profiles along the regions. All CPs appear as radial 
regions of steeper slope than ambient radial ranges. We de- 
termine their widths by fitting exponential functions to each 
profile inside and outside of the front as shown for the SE 
profile in Pig. [4] The steeper transition region between the 
radial ranges that are well- fitted is the cold front. We mark 
them by bold double lines in Pig. [4] The CP in the N is 
clearly detected in all three northern profiles. Towards the 
W, there is a double CP. The outer southern CP is more sub- 
tle, but clearly detectable in all three southern profiles. The 
CP positions and widths for each profile are summarised in 
Table |2] These widths are upper limits on the true front 
width, because they contain projection effects. IVikhlininl 
et al.| ( |2001| derived the X-ray brightness profile of an ellip- 
soidal gas cloud with a sharp outer edge. They showed that 
the brightness profile near the edge is approximately oc a/5, 
where d is the distance from the edge. The radial ranges we 
have identified as CPs and marked by double vertical lines 
in the profiles in Pig. ^ can be fitted by this function, sug- 
gesting that also the CPs in A496 are discontinuities in ICM 
density. 



Table 2. Positions and widths of cold fronts in A496 in profiles 
along regions marked in Fig. ^ 





distance of cold 


upper limit on 




front to center 


cold front width 


profile as labelled 


in kpc 


in kpc 


in Fig.|4] 


(in arcsec) 


(in arcsec) 


Nl 




7.6 (12) 


N2 


57 (90) 


7.6 (12) 


N3 




7.6 (12) 


W, outer 


57 (90) 


5(8) 


W, inner 


35 (56) 


5.7 (9) 


S 


158 (250) 


7.6 (12) 


sw 




13 (20) 


SE 




15 (24) 



2.7 Previous suggestions regarding the origin of 
the cold fronts in A496 

The observations of A496 have been interpreted differently 
by different authors. D07 argued that the spiral pattern 
of the cold fronts and their multiplicity are unlikely to be 
caused by a head-on merging remnant core and favoured 
the sloshing scenario. Comparing qualitatively to the sim- 



ulations of Ascasibar & Markevitch (2006), they derived 



subcluster orbit roughly from south to north. They spec- 
ulated that the subcluster found in the outskirts of A496 
towards the N-NW ( |Plin Krywult||2006| ) could have trig- 
gered the sloshing. 



Contrary, Tanaka et al. ( 2006 ) suggested that the non- 



uniform large-scale temperature and brightness residual dis- 



6 Roediger et al 



("-"offset) / "^Pc 




Figure 4. Width of cold fronts. Left: profile extraction regions perpendicular across the cold fronts. Right: Profiles along the regions 
shown on the left. The length coordinate along each profile is set off arbitrarily, it does not represent a distance from the cluster centre. 
In order to avoid confusion, we have multiplied/divided each profile as noted in the legend. The double vertical lines denote the position 
and width of each cold front, or more precisely an upper limit on its width. We determine the width by fitting exponential functions to 
each profile inside and outside of the front as shown for the SE profile. The steeper transition region between the radial ranges that are 
well- fitted is the cold front. We clearly identify the N cold front in profiles Nl, N2 and N3, and two fronts in the W profile. The fronts are 
more subtle in the three southern profiles, S, SE and SW, but clearly detectable. For comparison, the profile to the E does not contain 
a front. 



tribution is unlikely to be caused by minor- merger induced 
gas sloshing. Hence, they interpreted their southern CF as 
leading edge of a subcluster moving subsonically, and the 
northern one due to an oscillating central galaxy (see also 
Tittley Henriksen||2005t ). 

In this paper, we qualitatively and quantitatively verify 
gas sloshing as the origin of all observed features in A496 and 
infer an orbit similar to the prediction given in Sect. [23] on 
the basis of the position of the extended large-scale bright- 
ness excess in the N-NW. 



3 COMPARISON BETWEEN FIDUCIAL 
MODEL AND OBSERVATIONS 

In this section, we present the most likely merger scenario 
derived from comparing between simulations and observa- 
tions. We also discuss uncertainties and alternatives. 



3.1 The simulations 

For the sake of clarity, here we only briefly summarise our 
simulation method and fiducial simulations. A detailed ac- 
count of the simulation method, the initial cluster model, the 
subcluster model and orbit, and the synthetic observations 
are given in Appendix [B] 

We follow the approach of Rll and perform idealised 
simulations of a minor merger between a main cluster, 
here A496, and a smaller, gas-free subcluster. The ICM gas 



physics are described by the hydrodynamical equations. Ad- 
ditionally, the ICM is subject to the gravitational forces due 
to the main cluster and the subcluster. These are modelled 
by the adapted rigid potential approximation described and 



verified by [Roediger ZuHone ( 2011| ). Instead of modelling 
the DM distributions by the N-body method, we assume 
static potentials for both clusters, which speeds up the sim- 
ulations considerably. The subcluster potential is shifted 
through the main cluster along a test particle orbit. The 
simulations are run in the rest frame of the main cluster. As 
this is not an inert ial frame, the resulting inert ial accelera- 
tions are taken into account. 

The main cluster, here A496, is assumed to be initially 
in hydrostatic equilibrium and is set up according to ob- 
served ICM density and temperature profiles. Like A496, 
our model cluster is slightly elliptical and elongated along 
the grid ^/-axis. 



The subclusters are described as Hernquist halos ( Hern- 
quist||1990t . 



3.2 Summary of fiducial scenario 

Our comparison of simulation results and observations leads 
us to the following two fiducial runs that bracket the merger 
scenario: 

• Our fiducial subcluster has a mass of 4 X IO^^Mq and 
a scale radius of 100 kpc. 

• It crosses the main cluster along a diagonal orbit from 
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3.3.1 Spiral morphology and line-of-sight 

The spiral-like arrangement of the brightness excess and the 
CFs is a typical signature of gas sloshing seen along a line- 
of-sight (LOS) approximately perpendicular to the orbital 
plane of the subcluster. A LOS about parallel to the orbital 
plane would result in concentric arcs on opposite sides of the 



Figure 5. Fiducial subcluster orbits in the xy-plane. Like the sim- 
ulated maps in Fig.[l] we have rotated this plot by 15^^ clockwise, 
such that N is up and W is right. We mark the subcluster posi- 
tion in steps of 250 Myr as predicted from the test particle orbit 
and highlight its position at 0.75 and t = 1 Gyr for the "close" 
and "distant" subcluster orbit, i.e. with pericentre distances of 
100 kpc and 400 kpc, respectively. We note that the restricted 
potential approximation used here slightly over-predicts the cold 
front age and that the test-particle orbit neglects the subcluster's 
deceleration by dynamical friction, which is significant after peri- 
centre passage. Hence, we expect a current subcluster position 
considerably closer to the cluster centre, approximately 1 Mpc 
N-NW of the centre. We also mark the position of the subcluster 
found by |Fhn Krywultj ( |2006t . 



the SW to the N-NE, passing the cluster core SE of the 
centre at t = as depicted in Fig. |5] The two fiducial orbits 
have pericentre distances of 100 and 400 kpc, we refer to 
them as the "close" and the "distant" fiducial runs. 

• The simulation steps that match the observations best 
are 0.75 Gyr and 1 Gyr after pericentre passage for the 
"close" and the "distant" orbit, respectively. 

• Our line-of-sight is close to perpendicular to the orbital 
plane. 

In the following subsections, we give the reasons why 
we regard this configuration as the most likely. 



3.3 General features and merger geometry 

Along with the observed X-ray image, brightness residual 
maps and temperature map. Fig. [l] also displays the corre- 
sponding maps for our fiducial simulations, using the same 
colour scales. We have rotated our simulated maps by 15^^ 
clockwise to match the observed orientation of features. In 
this configuration, the fiducial runs reproduce all qualita- 
tive sloshing features simultaneously: the spiral pattern of 
the CFs and the brightness excess, the orientation of the 
brightness excess spiral, the large-scale brightness excess to- 
wards the N-NW, the strongest large-scale brightness deficit 
towards the E, and the appearance of the temperature map. 



cluster core (see Fig. CI), which is not the case here. 



3.3.2 Orientation on the sky 

The orientation of the subcluster orbit on the sky can be 
constrained by two properties: 

• the orientation of the brightness excess spiral, if its full 
extent is known, 

• and the orientation of the large-scale brightness asym- 
metry induced by sloshing. 

3.3.2.1 Fiducial orientation: In A496, the full extent 
of the brightness excess spiral is covered by the XMM- 
Newton observation. In the configuration described above, 
our fiducial simulation reproduces both properties simulta- 
neously. Thus, we derive a subcluster orbit from the SW 
to the N-NW with a pericentre SE of the cluster core, as 
sketched in Fig. |5] 

From their sloshing simulations for the Virgo cluster, 
Rll predicted the existence of a large-scale asymmetry as 
observed here. In Virgo, this asymmetry could not be con- 
firmed observationally, because the data did not extend far 
enough. Moreover, Virgo is still a very dynamic cluster and 
the sloshing-related asymmetry may never be disentangled 
from other perturbations. Here, in A496, we can for the first 
time confirm the connection between gas sloshing and the 
large-scale asymmetry, which is found in simulations and 
observations at the same position. 

3.3.2.2 Uncertainty of orbit orientation due to pre- 
ferred direction of sloshing in elliptical cluster: The 

elongation of an elliptical cluster introduces a preferred di- 
rection of sloshing along the axis of elongation, introducing 
an uncertainty in the orbit orientation. In order to assess 
this uncertainty, we have tested different orbit orientations 
in the xy-plane. The outcome of this test is demonstrated 
in Fig. |6] The morphology and positions of the innermost 
northern and southern CF are very similar to the observa- 
tion in both cases, reflecting the preference for the direction 
of sloshing. The large-scale residuals, however, differ. The 
morphology of the outer southern CF differs from the ob- 
servation in both cases, especially when the orbit is parallel 
to the y-axis. The brightness excess in the N-NW is similar 
to the observed one when the orbit is parallel to the x-axis, 
but then the observed E-W asymmetry is missing. The lat- 
ter is present when the orbit is parallel to the y-axis, but in 
this case the N-NW excess is much too weak. If we rotate the 
orbit even further clockwise such that the subcluster finally 
moves towards the N-NW, the outer large-scale brightness 
excess flips to the S, which disagrees with the observation. 
Hence, the large-scale structure clearly favours an orbit ori- 
entation comparable to our fiducial case. 

3.3.2.3 Higher age combined with orbit rotated by 
180"^: If we had data on only the central ^ 200 kpc, there 



8 Roediger et al 



IK. 



X-ray brightness/azimuthal average axis: z t=1000.14Myr 




M 



Figure 6. Large-scale brightness residual maps. In the top panel, 
the orbit is parallel to the x-axis, in the bottom panel parallel to 
the y-axis. The simulations are identical to the fiducial run oth- 
erwise. We have oriented the panels such that they match the 
observations. While the morphology and positions of the inner- 
most northern and southern cold front are very similar to the 
observation, the large-scale residuals are not. The morphology of 
the outer southern cold front differs from the observation in both 
cases, especially in the bottom panel. The brightness excess in the 
N-NW is similar in the top panel, but the observed E-W asym- 
metry is missing. The latter is present in the bottom panel, but 
the N-NW excess is much too weak. 



would be another degree of freedom. We could rotate our 
simulated maps by 180"^ and wait another couple of 100 
Myr. Then the major CF, now in the S, would have moved 
further south to replace the observed outer southern CF. 
The secondary, now northern CF would have moved out far 
enough to mimic the CF observed in the N, and a tertiary 
front in about the right position in the S would match the 
observed inner southern CF. Even though the overall mor- 
phology of these fronts would be somewhat more disturbed 
than observed, at t = 1.8 Gyr the secondary and tertiary 
front would simultaneously be at the correct radii. The ma- 
jor, now outer southern CF would fall 50 kpc short of the 
radius of the observed position. In this situation, we would 
not know the full extent of the central brightness excess spi- 
ral, nor the large-scale brightness excess, and thus would not 
be able to unambiguously infer the orbit orientation. How- 
ever, our knowledge of the brightness distribution outside 
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Figure 7. Cold front radii as a function of time. Cold front radii 
are derived from directional profiles in the sectors indicated in the 
top of each panel. The thin black horizontal lines mark the radial 
range of the observed cold fronts. The thin vertical line marks the 
inferred cold front age for the "distant" fiducial run, the dashed 
vertical line for the "close" fiducial run. We show the results of 
all simulation runs in light grey lines and highlight the fiducial 
"distant" run (black) and the fiducial "close" run (blue dashed). 
In the latter case, the fiducial subcluster passes the cluster core 
at a distance of only 100 kpc and considerably distorts the cool 
core. This is also the strongest impact case studied. Otherwise, 
the cold front radii depend mostly on time and only moderately 
on subcluster and orbit characteristics. The error bars indicate 
the width of the fronts in the synthetic observations. 



200 kpc clearly argues against this configuration and we can 
disregard it. 



3.4 CF radii and age 

With time, all CFs move outwards. This outward motion 
is mainly governed by the underlying potential of the main 
cluster and depends only moderately on the properties of 
the subcluster and its orbit. This is demonstrated in Fig.[7| 
where we plot the evolution of the CF radii towards different 
directions for several simulation runs. The same effect was 
already seen for the Virgo cluster (Rll). Thus, the cluster- 
centric radii of the CFs are a good tracer of their age, i.e. the 
time since the subcluster's pericentre passage. At a simula- 
tion time of 1 Gyr, our fiducial model "distant" simultane- 
ously reproduces the CF radii towards all directions within 
5 kpc. Given that our model uses a static cluster potential, 
this is a remarkably good agreement with the observations. 
The CF ages derived from different simulations vary between 
0.9 to 1.1 Gyr, and in each case the simulation simultane- 
ously reproduces the CF radii towards all directions with 
the accuracy stated above. The only exception is the case 
with the strongest impact, the fiducial "close" run, where 
the subcluster passes the cluster centre at a distance of only 
100 kpc and significantly distorts the cool core. Here the 
CFs move outwards somewhat faster, reaching the observed 
positions in the N, S, and W already at 0.75 Gyr. There is 
no CF towards the E. 

Our rigid potential method overestimates the CF age 
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by about 200 Myr (see Roediger ZuHone|2011 ), hence we 
conclude that the CFs in A496 are about 0.6 to 0.8 Gyr old. 
DOT estimated an age of about 0.5 Gyr from a qualitative 
comparison with the simulations of |Ascasibar Markevitchl 
(2006), which were aimed at a more massive cluster. 



X-ray brightness/azimuthal average axis: z t=1100.42Myr 



3.5 Contrasts across the cold fronts and 
subcluster mass/size/orbit 

Generally, the impact of more massive and/or more com- 
pact subclusters on orbits with smaller pericentre distances 
leads to stronger contrasts in all quantities across the CFs. 
Using this trend, Rll were able to constrain the subcluster 
mass within a factor of 2 for the case of the Virgo cluster. 
Attempting the same for A496, in Fig. |3]we compare direc- 
tional profiles for X-ray brightness, projected temperature 
and metallicity for simulations and observations. 

We plot results from three representative simulations: 
the two fiducial runs and a low-mass case with a subcluster 
mass of only IO^^Mq and a pericentre distance of 100 kpc, 
which should leave a milder imprint. As in the "close" fidu- 
cial run the CFs move out faster, we consider a timestep of 
0.75 Gyr for this case instead of 1 Gyr used in the other two 
cases. 

The two weaker impact cases achieve an excellent agree- 
ment with the observations within the inner 130 kpc. For the 
fiducial "close" run the temperature contrasts are slightly 
too large. As discussed in Sect. |2.5| the situation in A496 
is complicated by the fact that the CFs reside still inside 
the cool core of the cluster, where the CF jumps are in the 
same directions as the intrinsic gradients of the cluster, and 
thus blend with the general cluster profile. Consequently, 
the simulations show a strong degree of degeneracy. 

Only the regions outside the cool core towards the S and 
the E, i.e. in the region of the outermost southern CF, show 
distinguishable results. Here, subclusters with masses below 
2 X 10^^ Mq leave a negligible imprint and are incompatible 
with the observations. The best match in temperature and 
metallicity across the outer southern CF is achieved with 
the fiducial "close" run. In the fiducial "distant" case, the 
temperature and metallicity gradients are not quite strong 
enough, and the truth is presumably in between these cases. 
Thus, our comparison suggests a subcluster mass of 2 to 
4 X IO^^Mq. For comparison, the inner 150 to 200 kpc of 
A496 contain 4 x IO^^Mq, and the cluster mass within 1 
Mpc is 2 X IO^^Mq . A more precise constraint of the subclus- 
ter and orbit properties requires a full hydrodynamics+N- 
body treatment in the simulations and tighter observational 
constraints on temperature profiles and contrasts across the 
fronts. 



4 DISCUSSION 

4.1 Variations in initial cluster model 

4.1.1 Spherical models for A496 

Spherical cluster models for A496 reproduce the observed 
profiles for all quantities inside 150 kpc as well as the ellip- 
tical ones. Our estimate for the CF age marginally increases 
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Figure 8. X-ray brightness residual map for the initially spherical 
cluster. The sloshing spiral is more regular. The fiducial subclus- 
ter moved on a diagonal orbit (see black line). 



by 100 Myr at most. Outside 150 kpc, the spherical simula- 
tions slightly over-predict the X-ray brightness in E and W 
direction because they neglect the intrinsic ellipticity. 

The brightness excess spiral is, however, less elongated 
compared to the elliptical cluster and much more regular 
(Fig.|8|. With a diagonal SW-NE orbit orientation (see black 
line in Fig. [sj, these models also reproduce the large-scale 
brightness excess towards the N-NW, but the outer bright- 
ness deficit is located towards the SE. This difference to the 
elliptical models arises because the elliptical structure alone 
leads to a characteristic brightness enhancement towards N 
and S and deficits towards E and W. 

Overall, the elliptical cluster model results in a better 
match. 



4.1.2 Influence of initial temperature profiles 

From the observations, the temperature profile of A496 is 
not well-constrained. Hence we tested several possible fits to 



the observations (see Fig. Bl). Our results are not sensitive 
to the details of the temperature profile. The best match is 
achieved with profile 2, which is also the basis of our elliptical 
cluster model. 



4.2 The nature of the outer southern cold front 

In our simulations the outer southern CF never is a true 
discontinuity, while the observations suggest that it is. This 
is the only aspect the simulations do not reproduce. In 
the stronger merger simulation, the gradients across it are 
steeper in all quantities, but it is still continuous. Rll found 
the same structure in their Virgo sloshing simulations and 
named it cold fan to distinguish from a proper cold front. 
We verified that this is not an artefact of the rigid potential 
approximation employed here ( Roediger ZuHone|2011 ). A 
different orientation between the line-of-sight and the orbital 
plane would not make the feature appear sharper. 

Alternatively, we could alter our fiducial merger geom- 
etry and history in order achieve a true outer southern CF, 
but we would sacrifice the good match we achieve in all other 
properties: After about 4.5 Gyr after pericentre passage, the 
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now inner southern CF would have moved out far enough to 
reach the position of the observed outer front and could re- 
place it. At this stage, the northern CF would have reached 
a distance of about 270 kpc, where no CF is observed. The 
same is true for the western CF. An inclination of 77° be- 
tween the orbital plane and the line-of-sight would shorten 
the apparent distances of the CFs to the observed values, 
but only along one axis. At this strong inclination, the clear 
spiral pattern would be absent. The second option is to con- 
sider the stage at 1.8 Gyr after pericentre passage and rotate 
the merger geometry by 180"^, such that the previously outer 
northern front replaces the outer southern one, and the pre- 
viously inner southern one replaces the outer northern one. 
As discussed in Sect. |3.3.2.3| this configuration does neither 
match the overall brightness distribution, especially at large 
scales, nor all CF radii simultaneously. Thus, our fiducial 
configuration is the best match. 

The single minor merger scenario we suggest for the 
history of A496 is a very simple one, and an almost com- 
plete match between predicted and observed features can be 
regarded as a success. Very likely the history of A496 has 
been more complex, e.g. there may be remnant bulk flows 
from an earlier merger. The interaction of outer southern 
cold front/fan with this remnant perturbation may transfer 
the cold fan into a sharp cold front. 

4.3 Identifying the responsible subcluster 

Flin & Krywult f2006) analysed the galaxy distribution in 



A496 out to large radii and detected a subcluster or galaxy 
group towards the N-NW at the outer edge of their field-of- 
view. They applied their substructure detection method to a 
large number of clusters, and their general description of the 
method implies that each cluster is studied out to 1.5 Mpc 
from the cluster centre. Specifically for A496, they compare 
their result to the ones of Durret et al. (2000). This compar- 



ison reveals that in the case of A496 the field-of-view was 
larger, and that the subcluster in question is about 55 ar- 
cmin or 2.1 Mpc from the cluster centre towards the N-NW. 
We have marked the subcluster position in Fig.[5]along with 
the predicted position of our simulated subcluster, i.e. at 
1 Gyr after pericentre passage. 

Our prediction of the current subcluster position re- 
quires two more considerations. First, the rigid potential 
method slightly over-estimates the CF ages. Secondly, the 
employed subcluster orbit is a simple test mass orbit, which 
does not include dynamical friction. While this is accurate 
enough for the orbit prior to pericentre passage, dynami- 
cal friction will slow down the subcluster significantly after 
pericentre distance (see Roediger & ZuHone 2011). Both 



considerations require that we expect the subcluster at a 
much smaller distance to the cluster centre than naively ex- 
pected in our first estimate, i.e., around 1 Mpc towards the 
N-NE. While there is some intrinsic uncertainty in the di- 
rection of the orbit due to the ellipticity of the cluster, in 
Sect. |3.3.2.2| we have excluded orbits which bring the sub- 
cluster towards the N-NW. Furthermore, dynamical friction 
makes it unlikely that a subcluster that passed the cluster 
centre only about 1 Gyr ago or less can already have reached 
the observed distance of the Flin-Krywult subcluster. These 
reasons argue against this subcluster being responsible for 
the sloshing signatures in A496. 



The only other subcluster detected towards the N is 
LDCE0308 about 1.3 Mpc to the N-NW. While its distance 
is favourable, our simulations clearly favour a position to- 
wards N-NE instead of N-NW. Furthermore, this subcluster 
is detected in the 2MASS survey only and probably of low 
mass. 

This leaves the identity of the subcluster still an open 
question, because in the region where we expect the sub- 
cluster, there is no obvious galaxy concentration. However, 
it may be impossible to ever identify the responsible sub- 
cluster due to the effect of tidal forces. During the pericen- 
tre passage, a subcluster is tidally compressed. A few 100 
Myr afterwards, however, it suffers substantial mass loss due 
to tidal decompression or tidal stripping, and may well be 
dispersed and not recognisable as a compact structure any- 
more. 

Extrapolating from the simulations of [Ascasibar Sz\ 
[Markevitch| ( |2006| ), DOT speculated that the disturber should 
be rather massive with about 8 x IO^^Mq. Our simulations 
show that such a massive disturber very likely causes too 
strong an impact and destroys the cool core, but a less mas- 



sive on of 4 X 10^^ is sufficient. 



4.4 Origin of the disturbed inner cold fronts 

In high-resolution sloshing simulations without viscosity or 
magnetic fields, the sloshing cold fronts tend to be sub- 
ject to the Kelvin-Helmholtz instability (KHI) (e.g. Rll, 
ZuHone et al.||2011 ). The KHI can be suppressed by mag- 
netic fields aligned with the fronts ( ZuHone et al.|20lT ), and 
the fronts also appear smoother in a viscous ICM ( ZuHone 



et al. 2010). While in many clusters the CFs indeed ap- 
pear as smooth arcs free from signs of instabilities, in A496 
they show kink-like deviations from a smooth arc on length 
scales of about 20 kpc. In our "close" fiducial run, we re- 
produce KHI-induced distortions in the inner CFs that re- 
semble the observed ones, also the double CF towards the 
W. The length scale of the simulated distortions is even in- 
dependent of resolution (Fig. D3), indicating that the com- 
plex interplay of the curvature of the shearing layer and its 
outward motion set the instability length scale instead of 
the numerical viscosity. The presence of KHI puts an upper 
limit on the magnetic field strength in this cluster, because 
a mean tangential magnetic field, B, at the CFs stronger 
than fVikhlinin Markevitc"h||2002 ) 
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should suppress the KHI, where picM is the ICM pressure. 
Ma is the Mach number of the shear flow, and Thot and Tcoid 
are the temperatures at the warmer and colder side of the 



discontinuity, respectively. ZuHone et al. (2011 ) showed that 



gas sloshing typically amplifies the magnetic fields at the 
CFs by up to an order of magnitude, implying initial field 
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strengths of below a few /xG or magnetic pressure below 
several hundredths in A496, which is within observational 
limits for galaxy clusters in general. 

We note that the presence of substructure in CFs is not 
restricted to A496. Also the CFs in the centres of NGC 7618 
and UGC 12491, a pair of merging galaxy groups (Kraft 
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et al. 20061, show significant substructure in the form of 
kinks and wings (Kraft et al., in prep.). The same is true for 



the leading edges of M89 (NGC 4552, Machacek et al.|2006 ) 
and NGC 4472 (M49, |Kraft et al.|2011t , two Virgo ellipticals 
moving through the Virgo cluster. Thus, the presence of KHI 
at the CFs in some clusters and their absence in others may 
provide a sensitive probe for magnetic field strengths. 

Both, the presence of the instabilities in the inner CFs 
and the stronger contrast across the outer southern one 
favour the fiducial "close" run, but the cleaner match for 
the inner CF positions and contrasts across them favours the 
fiducial "distant" one. A distinction between both requires 
better observational constraints on the cluster potential and 
including the evolution of the DM components of both clus- 
ters in the simulation. Especially in the close encounter case, 
the central potential of A496 will deform somewhat, which 
will influence the exact positions of the CFs. 



4.5 Reconstructing the mass distribution in a 
sloshing cluster 



Markevitch et al. (2001) pointed out that the presence of 



the CF discontinuities in ICM density and temperature pro- 
flies can cause flawed results in cluster mass estimates. If 
the density and temperature proflles across a CF are used 
to derive the cluster mass proflle assuming hydrostatic equi- 
librium, the resulting mass proflle will have an unphysical 
discontinuity at the CF radius. 'Keshet et al. f2010) discussed 



that the tangential flows associated with sloshing naturally 
lead to centrifugal accelerations that add to the radial force 
balance. 

Here, we investigate to what extent the underlying clus- 
ter mass can be inferred from fully azimuthally averaged pro- 
flles. To this end, in Fig. [9] we plot the evolution of fully az- 
imuthally averaged (i.e. over 360"^) proflles for X-ray bright- 
ness, projected temperature and metallicity. 

Except for the central heating and brightness decrease, 
which are due to our neglect of radiative cooling, sloshing 
modifles the azimuthally averaged proflles by less than 10 
percent. The "close" flducial case, the distortions of the av- 
eraged proflles are weaker in the centre, but stronger at large 
radii. The inner cold fronts leave clearer traces in the "dis- 
tant" flducial run. Here, the modiflcation is systematic: in- 
side each CF, brightness and metallicity are enhanced and 
the temperature decreases, while outside each CF the op- 
posite is true. In order to detect this effect observationally, 
both, the brightness and the temperature proflle need to be 
known with at least 10 kpc resolution and with an accuracy 
of better than 10 percent. Given the scatter in the available 
data, the distortion of the proflles due to sloshing does not 
seem to introduce an observable difference and will not be 
a major error source for cluster mass proflles. The observed 
fully azimuthally averaged proflles for A496 show no trace 
of the CFs. 

Fig.[9]also shows that the accuracy of the averaged tem- 
perature proflle will depend on its derivation. Dividing the 




Figure 9. Evolution of fully azimuthally averaged profiles of 
X-ray brightness (top), projected temperature (third) and metal- 
licity (fifth), in our two fiducial models (see legend). The second, 
fourth and sixth panel show the relative difference to the initial 
state. The vertical double lines mark the cold fronts towards the 
N (dark blue) and S (light blue). For the temperature and metal- 
licity in the fiducial "distant" run, we show the results of two 
averaging approaches: either, we first derive a synthetic map and 
average azimuthally over annuli in this map (solid lines). Alter- 
natively, we calculate the integrals in Eqn. |B2] not only along the 
LOS for each pixel, but also over annuli, and perform the divi- 
sion in Eqn. |B2] afterwards (dashed lines). This corresponds to 
first binning the observed data over annuli into cumulative spec- 
tra and deriving the temperature and metallicity from them. The 
latter method leads to stronger imprints of the cold fronts in the 
profiles. The observed fully averaged profiles are also plotted here 
for comparison, but show no trace of the cold fronts. 
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X-ray data into annular rings and fitting the cumulative 
spectrum of each ring with a single temperature results in 
a stronger imprint of the CF in the profile (dotted coloured 
lines). Instead, profiles could be derived in different direc- 
tions and averaged azimuthally afterwards (solid coloured 
lines), which leads to a smaller deviation. Thus, for the pur- 
pose of mass reconstruction, the latter approach should be 
taken, while the former one will make CFs easier to detect 
already in averaged profiles. 



4.6 Metal redistribution 

In the bottom two panels of Fig. [9] we demonstrate that the 
metal redistribution introduced by gas sloshing is mainly 
an oscillation around the original metallicity profile. Hence, 
sloshing does not truly broaden the metallicity distribution. 
The same features were found in the Virgo cluster (Rll). 



5 SUMMARY 

We have presented a consistent scenario of minor-merger 
triggered gas sloshing in the slightly elliptical cluster of 
galaxies A496 by comparing observations to dedicated sim- 
ulations. We have deduced the following scenario as most 
the likely: about 0.6 to 0.8 Gyr ago, a subcluster of about 
4 X IO^^Mq passed between 100 and a few 100 kpc SE of 
the cluster core, moving from the SW to the N-NE. The 
subcluster should now be located about 1 Mpc N-NE of the 
cluster centre, but is probably already dispersed by tidal 
disruption and thus maybe impossible to identify. No obvi- 
ous subcluster is found in the galaxy distribution near the 
expected position. 

Our proposed scenario can explain almost all observed 
features of the ICM in A496: the spiral-like surface bright- 
ness excess wrapped around the cluster core and the asso- 
ciated asymmetries in temperature and metallicity maps, 
the positions of the cold fronts along the outer edge of the 
brightness excess spiral, the profiles for X-ray brightness, 
temperature and metallicity inside 150 kpc, and the large- 
scale distribution of brightness excess and deficit within the 
observed field of view, which extends out to 500 kpc. Only 
the outermost southern cold front is not a discontinuity in 
our simulations. 

For a small pericentre distance of 100 kpc, our simula- 
tions even reproduce the kinks observed in the inner cold 
fronts. This suggests that they are Kelvin-Helmholtz insta- 
bilities arising from shear flows along the cold fronts. This 
puts an upper limit of roughly 10 /iG or a magnetic pressure 
of 10% of the thermal pressure on the magnetic field parallel 
to the cold front. 

We have studied the effect of gas sloshing on fully az- 
imuthally averaged profiles for X-ray brightness, projected 
temperature and metallicity and found deviations from the 
original profiles to be below 10 percent. Thus, only in very 
deep data, traces of sloshing can be detected in fully aver- 
aged profiles, and only then sloshing could introduce errors 
in the mass reconstruction by violating the hydrostatic equi- 
librium. 
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APPENDIX A: COMPARISON OF 
OBSERVATIONAL DATASETS 

Detailed observations of A496 have been presented by sev- 
eral authors. Here, we use the datasets listed in Table [AT] 
In Fig. |Al| we compare directional profiles for X-ray bright- 
ness, projected temperature and metallicity. From the XMM 
archival data used by GIO and Lll, we extract spectra from 
concentric annuli to determine the projected temperature, 
metallicity and iron abundance profiles. First we fit the data 
with an absorbed APEC model in the 0.4-10 keV band with 
the hydrogen column density, uh , fixed at the galactic value 
to obtain temperature and metallicity (labelled "APEC"). 
Subsequently, the iron abundance has been determined us- 
ing a VAPEC model in the same energy band with uh first 
fixed (labelled "Fix") at the galactic value and then left free 
to vary (labelled "Free"). During the fitting procedure the 
abundances of O, Mg, Si, S, Ar, Ni were left free to vary. 
The ill-constrained hydrogen column density towards A496 
results in the systematic scatter in temperature and metal- 
licity (or iron abundance) between the different assumptions 
and also the difference to the GIO profiles. 

The X-ray brightness profiles from the D07 Chandra 
data fall significantly below the ones from the GIO XMM- 
Newton data outside ^ 70kpc, which is due to different 
methods of background subtraction. In the Chandra obser- 
vation, the whole field-of-view contained cluster emission, 
leading to an over-subtraction of the background. Therefore, 
for comparisons to simulations, we use the X-ray brightness 
profiles from GIO. The simulated temperature and metallic- 
ity profiles will be compared to the ones of GIO and our three 
recalculations, where we regard the range of observational 
results collected here as representative for their systematic 
error. 



APPENDIX B: SIMULATION METHOD AND 
INITIAL MODEL 

We follow the approach of Rll and perform idealised simu- 
lations of a minor merger between a main cluster, here A496, 
and a smaller, gas- free subcluster. The ICM gas physics are 
described by the hydrodynamical equations. Additionally, 
the ICM is subject to the gravitational forces due to the 
main cluster and the subcluster. These are modelled by the 
adapted rigid potential approximation described and veri- 
fied by Roediger & ZuHone|f2011 ). Instead of modelling the 
DM distributions by the N-body method, we assume static 
potentials for both clusters, which speeds up the simulations 
considerably. The simulations are run in the rest frame of the 
main cluster. As this is not an inertial frame, the resulting 
inertial accelerations are taken into account. We concentrate 
on minor mergers, where the subcluster is significantly less 
massive than the main cluster. A major merger would cause 
additional structure, e.g., destroy the cool core, which is not 
the case for A496. In A496, only in the inner ~ 15kpc the 
cooling time is shorter than 1 Gyr. Hence, in the outer re- 
gions of interest for our analysis, the cooling time is long 
and we neglect radiative cooling in our simulations. Finally, 
we only model the first core passage of the subcluster and 
the subsequent gas sloshing. 
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Table Al. Observational datasets for A496 used in this work. 



Reference 



instrument 



observing date 



net exposure time in ks 
(for XMM-Newton: 
(MOSl+MOS2+pn) 



Dupke &L White III (20031 



Dupke et al. {2oWf 
Tanaka et al. ( 2006 ) 
Lovisari et al.^ ( ^2011| ), Ghizzardi et al. (20101 XMM-Newton 



Chandra, ACIS-S3 
Chandra, ACIS-S3 
XMM-Newton 



10/2001 
07/2004 
02/2001 

08/2007, 02/2008 



8.7 

60 

14.3 + 14.0 + 7.4 

(59+61) + (59+60) 



(42+46) 



r / arcsec 
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r/ kpc 
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r / kpc 



Figure Al. Comparison of observed profiles of X-ray brightness (top panels), projected temperature (middle panels) and metallicity 
(bottom panels) towards N and S (left figure), and E and W (right figure), for different datasets, see legend and Appendix [A] For the 
orientation and range of the sectors, see Fig. |B4| The dashed vertical double lines mark the position of the cold fronts as stated in 
Table [U 



Bl Code 



B2.1 Spherical model 



Our simulations use the FLASH code (version 3.2, Dubey 
et aL||2009|. FLASH is a modular block-structured AMR 



code, parallelised using the Message Passing Interface (MPI) 
library. It solves the Riemann problem on a Cartesian grid 
using the Piecewise-Parabolic Method (PPM). The simula- 
tions are performed in 3D and all boundaries are reflecting. 
We use a large simulation grid of size 4x3x3 Mpc^ to ensure 
that no reflected waves reach the cluster centre during our 
simulation. In most simulations, we resolve the inner 32 kpc 
with Ax = 2 kpc and the inner 128 kpc with Ax = 4 kpc. 
The fiducial run uses a twice as good resolution. We per- 
formed resolution tests for several cases by rerunning with 
a twice and a four times as good resolution. In Appendix [P] 
we show that none of our results depend on resolution. 



B2 Model for A496 

The main cluster A496 is modelled either as a spherical or 
triaxial ellipsoidal cluster. 



Given the ICM density and temperature profile for A496 and 
assuming hydrostatic equilibrium, we calculate the gravita- 
tional acceleration due to the underlying main cluster po- 
tential as a function of radius. When the cluster is allowed 
to evolve for 2 Gyr in isolation, no sloshing or other modifi- 
cation of the ICM distribution is observed. 

The ICM density and temperature profile are fitted to 
observational data, which is summarised in Figure |B1| the 
first panel shows the radial electron density profile, ne(r), 
from the ACCEPT sample (ICavagnolo et al. 2009), from 



Kaastra et al.| ( |2004| (converted from proton density) , and 



the single-;^ fit given by Sivanandam et al. (2009). In the 



second panel, the grey-shaded regions display projected tem- 
perature profiles, Tproj (r), towards N, S, SE, W, E from GIO. 
The data points show azimuthally averaged projected tem- 



perature profiles from the ACCEPT sample and fr om [Kaas- 
tra et al.|(|2004t,|Snowden et al.| ( |200 8), Lll, and [SivanaiT 



dam et al. (2009). The third panel shows a collection of 
surface brightness profiles towards N, S, W, and E (DOT, 



GIO), and the azimuthally averaged profiles from ,Snowden| 
et aL|(|2008| and I Sivanandam et al.l (|2009|). 



We describe the density profile of our model cluster by 
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Figure Bl. Initial profiles of our model cluster and comparison to 
observations. The panels from top to bottom are electron density, 
temperature, are X-ray surface brightness. Black and grey solid 
lines show our deprojected model profiles, while dashed black and 
grey lines show projected quantities from our model. Coloured 
lines, shaded areas and data points are observational data from 
several works, see legends. Lines/symbols appearing in several 
panels are explained in the panel where they appear first. For 
details see Sect. Em 



a double- /3 profile. For the temperature, we utilise the ana- 
lytical formula 



Ti{r) 

step(x) 
drop(x) 



(m(r - ro) + To) x step(r ■ 
X drop(r — ro) with 

^-yi 

yi + 



■ro) 



(Bl) 



1 + exp [— (x ■ 



rs)/as 



which describes a linearly decreasing temperature in the 
outer part, a central decrease and an intermediate tem- 
perature enhancement. Given the scatter in the tempera- 
ture data we test three different temperature profiles for the 
model cluster. Our parameters for the density and temper- 
ature profiles are listed in Table [BT] 



Table Bl. ICM parameters for the cluster model. The fiducial 
model is based on temperature profile 2. 



Density: double /3 profile 

core radii ri,2/kpc: 

core densities poi,02/( gcm""^): 

/5i,2: 



Temperature: see Eqn. Bl 



m/(Kpc -*■): 
To/(107 K): 
ro/kpc: 
Vl 

Tsl kpc: 
as I kpc: 
r^/kpc: 
B 



profile 1 

5 


1.6 



33 
4.1 
2 



10 

5.64 • 10-26 
1 



profile 2 
-11.6 
4.7 


1.3 

12 

45 
29.6 
3.15 
2 



16 

6.6 • 10- 
0.42 



profile 3 



4.7 
-10 
1 



45 
2.5 
3 



B3 



Metal density: Hernquist profile, see Eqn. 

scale radius rn / kpc: 300 

core density pi:f /( gem- '^): 2.6 • 10-^^ 

metallicity fioor / solar: 0.33 



The black solid line in the top panel of Fig. |B1| displays 
our deprojected electron density profile. The thick dashed 
black line in the third panel displays the X-ray brightness 
profile of our model cluster, derived by projecting neA(T, z), 
where A(T, z) is the metallicity-dependent cooling function 
of |Sutherland & Dopita] ( |1993 ). Our density profile differs 
slightly from the fit given by Sivanandam et alT| ( |2009| ) , but 
results in a better fit to the X-ray surface brightness profiles. 
We experimented with the Sivanandam density profile, with 
both the original single-/^ version and with an additional 
central core. The resulting initial X-ray brightness profile 
runs along the upper range of the observed profiles shown 
in panel 3 of Fig. |B1| After sloshing, the X-ray profiles end 
up being too broad in all directions, hence we prefer our fit. 

The black and grey solid lines in the second panel of 
Fig. |B1 1 display our deprojected temperature profiles, while 
the black and grey dashed lines present the corresponding 
projected temperature profiles. These are derived following 



Mazzotta et al. ( 2004 ) by calculating a weighted average of 



the temperature along each LOS using the weights 



W = n IT 



3/4 



(B2) 



The top panel of Fig. |B2] demonstrates the observed, az- 
imuthally averaged metallicity in A496 from Snowd en et al.| 
(|2008 ), Lll (iron abundance instead of metallicity), and 



Sivanandam et al. (2009). Together with the ICM density 



profile, this translates into a metal density profile (bottom 
panel), which we fit with a Hernquist profile ( |Hernquist| 
p9Q. ), 



PFe(r) 



27r r (r/rH + lf 



(B3) 



Furthermore, we impose a metallicity floor at 0.33 solar. 
The resulting metal density profile is shown in the bottom 
panel of Fig. |B2| The parameters for the core density, pH, 
and scale radius, rn, are given in Table [BT] The resulting 
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Figure B2. Radial profiles of metallicity (top) and metal density 
(bottom): we show our projected and deprojected model and ob- 
servational data, see legends. Correctly speaking, the data from 
[Lovisari et al.| ( [2011| ) display the iron abundance instead of metal- 
licity, leading to the differences near the centre. 



deprojected and projected, emission-weighted metallicity is 
shown in the top panel in comparison to observational data. 



B2.2 Elliptical model 

A496 is sUghtly elUptical along the N-NW to S-SE direction 
(TOG, Lagan a et al. 201QJ . This feature can also be derived 
from a close inspection of the observed directional surface 
brightness profiles in the third panel of Fig. |Bl| The average 
brightness along the E-W axis declines somewhat faster than 
the one along the N-S axis, which indicates an elongation 
along the N-S axis. 

The spherical cluster described in Sect. |B2.1| can be 
made triaxial while maintaining hydrostatic equilibrium by 
replacing the spherical radius, r ■ 
elliptical radius. 



Y^x^ + + z^, with the 



p2 y2 £1 



&2 



(B4) 



in the double-/^ profile for the ICM density and the temper- 
ature profile in Eqn. |B1[ We use a = 1.2 and h — c — \^ 
which makes the cluster elongated along the ^-axis. For our 
elliptical cluster we use temperature profile 2. In Fig. [2] we 
compare the observed and modelled surface brightness pro- 
files for separate directions. 



B3 Subcluster and orbits 

We test several combinations of subcluster masses, scale 
radii and orbits. Once we know the main cluster's poten- 
tial, we calculate the orbit of a test mass moving through 
this potential. During the course of the simulation, the sub- 
cluster potential is shifted along its orbit through the main 
cluster. The time normalisation is chosen such that pericen- 
tre passage happens at t = 0. At the start of each simulation. 



Table B2. List of simulation runs, stating the subcluster mass, 
scale radius, pericentre and apocentre of its orbit. 

^sub/lO^^M0 dsub/kpC fimin/kpC cimax/ MpC 
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10 



the subcluster is placed at its orbit 1 Gyr prior to pericen- 
tre distance. We stop the simulations at t ~ 1.5 Gyr after 
core passage, shortly after a good match with the observed 
CFs is reached. This moment is well before the second core 
passage of the subcluster. 

The gravitational potential of our subclusters is de- 
scribed by a Hernquist halo ( Hernquist]|1990 ) . We vary its 



Mq and be- 



mass and scale radius between 0.5 and 4 x 10 
tween 50 and 200 kpc, respectively (see Table B2). A mass 
of 4 X IO^^Mq equals the mass of the inner 150 to 200 kpc 
of A496. In our cluster model, A496 contains 2 x IO^^Mq 
within 1 Mpc, thus our simulations are for mass ratios above 
5. At the upper end of our mass range, we are restricted by 
the condition that we want to model minor mergers only. 
The lower end of the mass range is determined such that 
the subcluster still leaves an imprint that is comparable to 
the observations. 

We assume our subclusters to be gas- free. Rll have 
shown that in the Virgo cluster a subcluster sufficiently 
massive to cause the observed CFs can be completely ram- 
pressure stripped before pericentre passage. The ram pres- 
sure in A496 in general exceeds the one in Virgo by at least 
a factor of 2, and we study the same subcluster mass range 
here. Thus, also A496 should be able to ram-pressure strip 
the subclusters under consideration. Also the observations 
do not show any sign of gas ram pressure stripped from a 
companion subcluster. 

We test pericentre distances of the subcluster orbit of 
100 and 400 kpc, and apocentre distances of 3 and 10 Mpc, 
leading to different velocities during pericentre passage. Fur- 
thermore, we vary the orientation of the orbits as described 
in Sect. 13.3.2.21 



B4 Synthetic observations 

From our simulations, we calculate synthetic maps and ra- 
dial profiles. For the X-ray brightness, we project n^A(T, z) 
along the LOS, where A(T, z) is the metallicity-dependent 



cooling function of Sutherland & Dopita ( 1993J . 

For each X-ray image, we calculate a brightness residual 
map by dividing the image by its azimuthal average, thus 
highlighting deviations from circular symmetry. This proce- 
dure requires the identification of the cluster centre around 
which the azimuthal average is taken. The details of the 
resulting residual map, especially near the cluster centre, 
depend somewhat on the positioning of the cluster centre. 
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Figure B3. Dependence of brightness residual maps on centre 
determination. The diamond marks the centre of the BCG in 
A496, the cross marks the centre used for the residual map. It is 
offset by 5.6 arcsec along the diagonals. The colour scale is the 
same as in Fig. [l] 




X-ray brightness/azimuthal average axis: z t=1000.23Myr 




Figure CI. Brightness residual maps along different LOSs: syn- 
thetic images for fiducial run at t = 1 Gyr along three grid axes. 
The elliptical cluster is elongated along the y-axis, the subcluster 
orbit is in the xy-plane, approximately from SW over SE to N-NE. 
Top left: along the 2;-axis, which is perpendicular to the orbital 
plane. Bottom left: along the y-axis. Top right: along the x-axis. 
Bottom right: residual map derived from the Chandra image of 
DOT. 

APPENDIX C: OTHER LINES-OF-SIGHT 

In Fig. |Cl| we show brightness residual maps along different 
LOSs. The CPs can be detected in all cases. The spiral- 
like arrangement of brightness residuals and CPs is only 
seen along the LOS perpendicular to the subcluster orbit. 
Por LOSs parallel to the orbital plane, CPs and brightness 
residuals take the appearance of staggered arcs. Hence, we 
conclude that in A496 we are looking approximately per- 
pendicular onto the orbital plane of the subcluster. 



Figure B4. Directional profiles are averaged over the sectors 
marked in this graph. See legend for different datasets. For 
our work, we consistently use the sectors marked with "sector 
N/W/S/E". 



as shown in Pig. |B3| However, the overall structure of the 
residual map including the spiral-shaped brightness excess, 
the northern CP and the brightness excess towards the S is 
detected robustly. 

Por the projected temperature, we employ the weight- 
ing scheme of Mazzotta et al. (2004) as explained above 
(Eqn. B2). The projected metallicity is the emission- 



weighted metallicity averaged along each LOS. 

Prom these synthetic maps, we calculate directional pro- 
files which are averaged over a 30° azimuthal range. We 
use the same azimuthal sectors as in GIO and as shown in 
Pig.|B4l 



APPENDIX D: RESOLUTION TEST 

We have performed resolution tests for several runs by vary- 
ing the resolution by a factor of 4. In Pig. |D1| we show ex- 
amples of brightness residual maps for our the run with the 
fiducial subcluster on an orbit parallel to the ^/-axis with 
pericentre distance 400 kpc, for three different resolutions. 
In Pig. |D2| we show projected temperature profiles for the 
same case at two different epochs. Only at late times and 
only in the inner 30 kpc, these profiles differ slightly for the 
low-resolution case. In Pig. |D3| we show X-ray and brightness 
residual maps for the fiducial "close" run. The distortions of 
the inner CPs are largely independent of resolution. 
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Figure Dl. Resolution test: brightness residual maps for our fiducial subcluster on orbit parallel to y-axis at t = 1 Gyr. The resolution 
increases by a factor of 2 from left to right between neighbouring images. 
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Figure D2. Resolution test for same case as in Fig. |Dl| projected 
directional temperature profiles for three different resolutions at 
0.8 Gyr and 1.3 Gyr after the pericentre passage. The resolution 
increases by a factor 2 from LR to HR and again to HHR. 



This paper has been typeset from a TJ]X/ I^T^X file prepared 
by the author. 
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Figure D3. Resolution test: brightness residual maps for our fiducial run "close" with pericentre passage at 100 kpc. The resolution 
increases by a factor of 2 from left to right between neighbouring images. In the rhs image, the northern CF is resolved with 1 kpc. The 
distortions of the inner CFs are largely independent of resolution. 



